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Relationships  between  Volume,  Surface  and  Line  Distributions  of  Vorticity, 
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Summary 

v 

Several  mathematical  theorems  are  derived  which  demonstrate  the  equivalence  of 
continuous  volume  distributions  of  doublicity,  vorticity  and  source,  and  show 
furthermore  that  their  influence  may  be  expressed  purely  in  terms  of  continuous 
surface  distributions  of  these  quantities  over  the  closed  boundary  of  the 
volume . 


These  general  theorems  may  then  be  particularised  to  'sheets'  of  singularities 
distributed  over  non-clos ed  surfaces;  amongst  a  number  of  examples,  the 
special  cases  of  the  velocity  fields  induced  by  source,  vortex  and  doublet 
sheets  are  considered,  which  under  certain  circumstances  are  equivalent  to  each 
other  and  reduce  to  simple  line  integrals. 

These  theorems  are  expected  to  have  some  application  in  aerodynamic  problems 
involving  the  interaction  between  irrotational  incompressible  flow  regions  and 
regions  of  more  general  flow  such  as  those  arising  in  aerodynamic  wakes  and  in 
the  jet-in-crossflow  problem,  and  to  be  of  assistance  in  the  development  of 
improved  surface  singularity  techniques.  -  .  ^  , 
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INTRODUCTION 


It  is  well  known  that  if  the  potential  in  a  region  obeys  Laplace's  equation 
everywhere  within  that  region,  then  that  potential  can  be  expressed  entirely 
in  terms  of  surface  integrals  involving  the  value  of  the  potential  and  its 
normal  derivative  everywhere  on  the  closed  boundaries  of  that  region.  This 
fact  is  exploited  in  practical  aerodynamic  calculation  methods  ('panel 
methods')  which  allow  the  potential  and  the  velocity  everywhere  on  and 
outside  the  surface  of  a  configuration  to  be  determined  nominally  exactly 
by  evaluating  one  or  the  other  of  these  functions  at  the  configuration 
surface,  such  that  the  appropriate  boundary  conditions  are  satisfied. 

However,  in  cases  where  the  governing  equations  in  the  region  of  interest  do 
not  reduce  to  Laplace's  equation  (such  as  flows  which  are  compressible  or 
rotational)  the  expression  for  the  velocity  at  any  point  in  such  a  region 
involves  not  only  surface  integrals  but  also  volume  integrals  which  in 
general  extend  over  the  entire  non-Laplacian  volume.  Consequently  schemes 
treating  such  cases  involve  a  significantly  greater  amount  of  computation  and 
currently  are  restricted  to  geometries  much  simpler  than  those  which  'panel 
methods'  are  able  to  handle. 

Problems  in  which  it  is  required  to  compute  the  flow  inside  regions  where 
volume  sources  or  vorticity  may  be  considered  to  be  present  (e.g.  compressible 
or  rotational  flow)  will  not  be  considered  here.  In  many  problems  of  interest, 
however,  such  regions  may  be  embedded  in  a  flow  which  is  otherwise  source-  and 
vorticity-free;  examples  include  the  flow  outside  a  thick  wake  (and  its 
'rolled  up'  core  of  rotational  fluid)  behind  a  lifting  wing,  or  the  flow  out¬ 
side  a  jet  issuing  at  some  angle  into  an  otherwise  irrotational  stream. 

In  such  cases  it  is  often  required  (sometimes  as  part  of  an  iterative  scheme) 
to  compute  the  effect  produced  in  the  outer  (Laplacian)  region  by  the  sources 
and/or  vorticity  used  to  model  these  embedded  regions.  It  is  shown  herein 
that  this  external  effect  may  instead  be  computed  from  surface  distributions 
over  the  boundary  of  these  embedded  regions,  and  that  in  certain  special 
cases,  when  the  embedded  region  may  be  modelled  as  an  infinitesimally  thin 
volume  (i.e.  a  sheets  line  distributions  may  be  used.  The  theorems  allow  the 
singularities  (vortices,  sources,  doublets)  to  be  chosen  which  are  computationally 
the  most  convenient  for  the  problem  in  hand. 


j  Accession  For 

NT IS  GRA&I 

jt 

DTIC  TAB 

Unannounced 

□ 

Just  if  teat  ion _ 

By _ _ _ 

j  Dlstrll:  it’  u>/ 

i  Availntij.it,*  C 

t.  V  «"l  i  t  ririi 

jDlst  1  Sp  i-jc  iai 

\&±\. _ 


-3- 


2. 


BASIC  CONVENTIONS  AND  INFLUENCE  EXPRESSIONS 


In  the  following,  a  general  vector  J  denotes  the  quantity  ( Ax  t 

where  t,  J,  1c  are  the  unit  vectors  of  the  cartesian  axes  x,  y,  z.  In 

particular  the  vector  r  denotes  the  vector  drawn  from  a  point  Q  to  a  point  P. 

i.e.  <r  =  +  Op  -  *■ 

or  T*  =  Cxp  +  Cz>.  -**)x 

The  operators  curl  a  (or  ^x) ,  divft  (or  V,*)  and  grada  (or  )  denote  the  usual 
vector  differential  operators,  the  differentations  being  defined  with  respect 
to  the  coordinates  of  the  point  Q.  The  operator  gradp  (or  )  implies 
differentiation  with  respect  to  the  coordinates  of  P. 


It  can  easily  be  seen  in  particular  that: 

(1)  $radp  (  -jr)  =  -frrad^C-r)  =  “ 

Use  will  be  made  of  the  standard  vector  identities: 


(2) 

div  (sA) 

s  A  .  grad,  s  +•  s  dtv  A 

(2a) 

cwrL  (sA) 

5  (^rad-s)xA  +  S  curt/) 

and 

(3) 

(B*C)x  D 

=  (6.5)  c  -  (c.5)  I 

where  s  denotes  a  distributed  scalar  function,  A  a  distributed  vector 
function  and  B,  C  and  D  are  general  vectors;  Gauss’  divergence  theorem  will 
also  be  used: 

( k )  JJf  div^  A  dJl  =  JJ  A  .  Tt  d  S 

ji  s 

where  A  is  a  continuous  vector  function  defined  at  every  point  Q  within  the 
arbitrary  volume  ft  and  n  denotes  the  outward  normal  to  the  surface  S  bounding  fL  . 


Consider  the  effects  induced  at  a  point  P  by  an  elemental  singularity  (source, 
doublet  or  elemental  vortex  filament)  located  at  a  point  Q;  this  may  be  an 
element  of  a  continuous  line,  surface  or  volume  distribution. 


The  elemental  potential  d$s  and  velocity  dVs  induced  at  P  by  an  elemental 
source  of  strength  do-  located  at  Q  Eire  given  by: 

(5)  d$j  =  do'  and  dV,  =  gradp  C d$s)  =  dcr 

The  elemental  potential  3<pd  and  velocity  induced  at  P  by  an  elemental 
doublet  dp  located  at  Q  (with  dp.  -  a  d/u.  ,  where  a  is  a  unit  vector  along  the 
doublet  axis,  the  positive  direction  of  the  axis  being  from  the  'negative' 
to  the  'positive'  end  of  the  doublet)  are  given  by: 


(6) 


VK  p-adp  (.+-)'  dp  =  -fa 
fa  fradr(dfa)  =  ~fa  f™*'?  ( 
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The  elemental  velocity  dVv  induced  at  P  by  an  elemental  vortex  filament  d  P 
located  at  Q  (the  positive  direction  of  the  vector  dT  being  associated  with  a 
clockwise  rotation  about  the  filament)  is  given  by: 

.  —  _  _j_  r  x  dr 

( 7 )  a  -f-3 

It  may  be  noted  by  comparing  (5)  and  (7)  that  the  velocity  field  induced  by  a 
unidirectional  Y°rt_icity  distribution  (line,  surface  or  volume )  may  be  derived 
from  that  induced  by  a  source_ distribution  with  the  same  spatial  density 
variation,  simply  by  computing  the  vector  cross-product  of  the  source- induced 
velocity  with  the  unit  vector  defining  the  vorticity  direction.  The  effect 
of  vorticity  distributions  of  variable  direction  may  simply  be  constructed  by 
compounding  three  different  unidirectional  vorticity  distributions  (even 
though  each  of  these  in  isolation  may  be  physically  impossible).  This 
subject  is  discussed  in  depth  in  ref.  1. 
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3. 


INTEGRATION  OF  ELEMENTARY  CONTRIBUTIONS 


3.1  Distribution  of  Vorticity  on  a  Closed  Surface 

3.1.1  Velocity  Induced  at  an  External  Point 


Consider  the  velocity  Vp  induced  at  a  point  P  due  to  a  continuous  surface 
vorticity  distribution  of  variable  density  T,  lying  everywhere  on  the  closed 
surface  S  of  an  arbitrary  volume  il  (see  Fig.  1  (a));  P  lies  outside  this 
volume.  Since  Tj  at  any  point  of  the  surface  must  be  perpendicular  to  the 

local  surface  unit  outward  normal  n,  we  can  write  Ts  -  ,  where  yd  is  some 

variable  vector  field  function  defined  at  the  surface  S. 


Using  (7)  the  velocity  induced  at  P  is  given  by: 

(8)  471  Vp  =  jj(R  *  n)  x  b  dS 

s 

where  b  is  an  abbreviation  for  the  variable  vector  r/r3  (i.e.  -  gradp(i=)). 
Now,  using  (3): 

(yl  »  n)  k  b  =  (jH.'S)'fx  —  (-n  .  5  )y<Z . 

Inserting  this  in  (8)  and  considering  only  the  component  in  the  direction  x: 


4TiVp.t  =  Jfj f  [(/!.■£  )nx  -  (£.‘S)d,JdS 

s  ' 

=  // [*  (yd.E)]’n  <JS  -  J/CyC^F).  n  dS 
s  s 

which  becomes,  using  (4): 

(9)  4TT  Vp  .  £  =  fff  c)iv9  [i(/Z.'5)]dn  -  ffj 

■G.  JX 

the  integration  extending  over  all  points  Q  within  the  volume  ft. 

Now  using  (2)  : 

[Z  (yU.F)j  =  i  .  +  ( jl .  I  )  odv«  £ 

=  t  (yiZ.^)  etiv,  1  s  O 


and 

[yd*F] 


£ .  yd,  +  /K  h 
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since 


AlSO:  b  .©HKlaii*  =  b*  ^21  +  k  +  U  3^- 

Inserting  these  last  two  equations  into  (10)  and  rearranging  gives: 

*  ♦*<**-&) 
which  by  inspection  is  equivalent  to: 

(12)  4TTVp.t  =  [prad(f(/l.'5}]ciA.  +JJJ4 .  [E  xcturl^ji]  dA. 

If  the  process  following  equation  (8)  is  repeated  for  the  components  in  the 
y  and  z  directions,  it  follows  that  we  finally  obtain  from  (12): 

(13)  4TC  VP  =  fff  -  ^dF(^S)  d/l  -  JfJ-s  veuA^p.  AA  . 

-a  12 
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—  _  *3 

By  comparing  with  equations  (6)  and  (7)  and  remembering  that  b  =  r/r  it 
can  be  seen  that: 

kIJJ-r*L'Qc-V‘ul 

is  the  velocity  induced  at  P  by  a  volume  doublet  distribution  throughout  the 
volume  Jl  ,  of  variable  density  p  ,  (see  Fig.  1  ( b ) ) ,  and  that 

*  ^V1  dil 

is  the  velocity  induced  at  P  by  a  volume  vorticity  distribution  of 
variable  density  %  =  curl^,  JL  .  Since  the  equations  are  derived  for  an 
arbitrary  volume  doublet  distribution  p  and  an  arbitrary  external  point  P,  it 
follows  that  these  relationships  between  the  equivalent  doublet  and  vorticity 
densities  apply  in  a  local  sense  for  any  point  Q  within  Jl  or  on  S,  and  not 
only  in  an  integral  sense. 

The  volume  vorticity  distribution  Tv  is  thus  everywhere  solenoidal  since 
div^  Yv  =  div^  carl ^  ft  =  0  (i.e.  vorticity  is  neither  created  nor  destroyed 
at  any  point  within  11).  It  can  furthermore  easily  be  shown  by  constructing 
local  cartesian  system  [  %  ,->2  ,n)  at  any  point  on  S  that  the  vorticity  flux 
-n.curl p  dfd entering  an  elemental  area  dS  =  d\dr^  from  the  interior  of  Jl  is 
exactly  equal  to  the  flux  of  the  surface  vorticity  leaving  the  edges  of  dS  . 

It  follows  that  the  surface  vorticity  px'n  is  physically  meaningful  in  isolation 
only  in  the  special  case  cur  Ip  =  0;  unless  this  is  the  case,  the  surface 
vorticity  cannot  be  replaced  by  surface  doublicity  in  the  manner  discussed  later. 

By  equating  the  right-hand  sides  of  equations  (8)  and  (13)  and  rearranging  the 
terms,  we  obtain  the  following  theorem: 

Theorem  1  "The  velocity  induced,  by  a  volume  distribution  of  doublets 
of  arbitrary  volume  density p  ,  at  a  point  lying  outside  that  volume3 
is  identical  to  the  velocity  induced  at  that  point  bu  a  volume 
distribution  of  vorticity  of  volume  density  K  =  carlp  throughout  that 
same  volume ,  together  with  a  surface  vorticity  distribution  7^  on  the 
surface  of  that  volume ,  whose  surface  density  and  direction  are  given 
by  Ts  where  pt  denotes  the  local  value  of  /t  infinitesimally 

inside  that  surface3  and  n  denotes  the  local  outward  unit  normal  to 
that  surface". 

It  may  be  noted  that,  if  the  volume  doublet  distribution  is  known,  the 
velocity  induced  at  an  external  point  may  be  evaluated  explicitly  either 
directly  in  terms  of  the  volume  doublet  distribution  or  in  terms  of  the 
equivalent,  volume-plus-surface  vorticity  distribution.  In  the  special 
case  where  the  volume  vorticity  has  a  density  of  zero  (i.e.  curl/i  =0), 
this  velocity  may  be  obtained  purely  from  a  surface  integral.  This  statement 
requires  modifying  when  the  point  at  which  the  velocity  is  required  lies 
inside  the  volume  il . 


-8- 


3-1.2  Velocity  Induced  at  an  Internal  Point 

In  the  case  where  the  point  P  lies  within  the  volume  JI  the  value  of  4=  approaches 
infinity  for  points  Q  in  the  immediate  vicinity  of  P;  at  P,  is  no 

longer  defined,  and  the  steps  leading  from  equation  (9)  to  equation  (10)  which 
assume  that  V<j  (4=)  =  0,  are  no  longer  valid.  This  problem  may  be 
circumvented  by  constructing  a  small  sphere  of  radius  g,  about  P  (its  surface 
being  denoted  by  Z  and  its  volume  by  or  )  and  excluding  this  sphere  from  the 
region  of  integration;  the  point  P  may  then  be  considered  to  lie  outside 
the  volume  (/1-v),  and  the  above  equations  are  then  still  applicable  for  this 
volume  and  the  associated  surface  (S  +2).  The  equations  relating  to  the  complete 
volume  Jl  may  then  be  obtained  by  augmenting  the  right-hand  sides  of  equations 
(10)  and  (12)  by  the  expression 

[ftp*  dJl] 

which  was  zero  for  the  case  where  P  lay  outside  SI  .  Now  as  £-»0  the  function 
becomes  progressively  more  uniform  within  v  and  in  the  limit  the  above  expression 
may  be  written  ~fi^  $  ^LV<t  f>  dSi  .  By  Gauss'  theorem  (equation  (h))  this 
then  becomes  -/*■■<  JTE.  -d'  dS.  with  it' the  outward  unit  normal  to  the 

surface  £.  Since  r  is  defined  positive_from  Q  to  P  (i.e.  pointing  towards 

the  centre  of  this  sphere),  the  quantity  b.n'  is  equal  to  -\/tL  for  all  points 
Q  on  Z  so  that  the  above  surface  integral  reduces  to  hUp*  .  Consequently 
equation  (13)  becomes  for  a  point  P  inside  Jl  : 

(lb)  4irYp  =  fff-  f^ad.  C/Z.'B)  All  -  fjj-t  *  curl^p.  dJl  +  i+n/z  . 

Jl  Jt 

It  can  thus  be  seen  that  for  the  case  where  the  point  P  lies  inside  the 
volume  Si  ,  the  following  theorem  is  obtained: 

Theorem  la  "The  velocity  induced  by  a  volume  distribution  of  doublets 
of  arbitrary  volume  density  p  ,  at  a  point  which  lies  within  that  volume , 
is  identical  to  the  velocity  induced  at  that  point  by  a  volume 
distribution  of  vorticity  (of  volume  density  %-  curl p  throughout  that 
volume ),  together  with  a  surface  vorticity  distribution  on  the 
surface  of  that  volume ,  (whose  surface  density  and  direction  are  given 
by  %  where  /Zt  denotes  the  value  of  jui  infinitesimally  inside 

that  surface ,  and  ri  denotes  the  outward  unit  normal  to  that  surface ) , 
plus  a  term  which  is  locally  equal  to  the  vector  value  -Ju". 

It  may  be  noted  that  if  the  volume  distribution  p  is  known,  then  the 

velocity  may  be  evaluated  either  directly  in  terms  of  the  volume  doublet 

distribution  or  in  terms  of  the  equivalent  volume-plus-surface  vorticity 

distribution,  augmented  by  the  local  term  -Ju  .  In  the  special  case  where 

the  volume  vorticity  has  a  density  of  zero  (i.e.  curl^u  =0),  this 

velocity  may  be  obtained  purely  from  a  surface  integral  plus  the  local  term  -/T. 


3.2  Source  Distribution  on  a  Closed  Surface 

Referring  again  to  Fig.  la,  consider  the  potential  <pp  induced  at  a  point  P 
due  to  a  continuous  surface  source  distribution  of  variable  density  oj  on 
the  closed  surface  S.  Suppose  that  the  scalar  function  o,  is  expressed 
as  oj  =/Z-ri  where  p.  is  again  some  vector  field  function  defined  at  the 
surface  S. 
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Using  (5)  the  potential  induced  at  P  is  given  by: 

(15)  =  -JJ^-  iS 

s 

which,  using  (1)  may  be  expressed: 

4-Tt  <|>p  =  ~  fff  (£?)  6A 

SL 

and  which,  using  (2),  further  reduces  to: 

UK$r  =  ~[\\fL-Trad-<i  &)*&  "  J7T  dJI 

jr 

or,  using  (l),  to: 

(16)  4TC  <pr  =  JJJ  (4)  •/!  -  f/r  Jiv  IZdil. 

A  Jl 

By  comparing  with  equations  (6)  and  (5)  it  can  be  seen  that 

ht  fff  v™*?  Ci)'A  <*a 

ji 

is  the  potential  induced  at  P  by  a  volume  doublet  distribution  of  variable 
volume  density p.  throughout  Jl  and  that 


is  the  potential  induced  at  P  by  a  volume  source  distribution  of  variable 
volume  density  <yv  =  -div^u  throughout  JL . 

It  can  be  seen  that  none  of  the  above  steps  precludes  the  case  where  P 
lies  inside  A  ,  i.e.  equation  (l6)  is  equally  valid  for  points  lying  inside 
or  outside  Jl  . 

Once  again,  the  equations  are  valid  for  an  arbitrary  volume  doublet 
distribution  and  for  any  point  P;  consequently  the  relationships  between 
the  equivalent  source  and  doublet  distributions  apply  in  a  local  sense  for 
any  point  Q  in  the  respective  distributions,  and  not  only  in  an  integral 
sense. 

The  volume  source  distribution  is  everywhere  irrotational  since 
curl^o'y  =  curl^  div^^I  =  0  (i.e.  vorticity  is  neither  created  nor  destroyed 
at  any  point  within  A)-  Furthermore  it  follows  immediately  from  Gauss' 
theorem  (equation  (L))  that  the  total  volume  integral  of  the  source 
density  -  divp  is  equal  and  opposite  to  the  total  surface  integral  of  the 
source  density  /Z.n  . 

By  applying  the  operator  gradp  to  both  sides  of  equations  (15)  and  (16), 
the  corresponding  equations  for  the  velocity  Vp  are  obtained.  By  equating 
the  right-hand  sides  of  these  equations  and  real  ranging  the  terms,  the 
following  theorem  is  obtained: 

Theorem  2  "The  potential  and  velocity  induced  at  any  point  P  hy  a 
volume  distribution  of  doublets  of  arbitrary  volume 
density  ju.  are  identical  to  those  induced  at  the  same 
point  by  a  volume  distribution  of  sources  of  volume 
density  oj  *  -  divp  throughout  that  same  volume ,  together 
with  a  surface  source  distribution  eg  =  /v-n  on  the  surface 
of  that  volume ,  where  Ju. €  denotes  the  local  value  of  p 
infinitesimally  inside  that  surface ,  and  -n  denotes  the  local 
outward  unit  normal  to  that  surface". 


Since  Theorem  1  (and  la)  and  Theorem  2  are  stated  for  arbitrary  volume 
doublet  distributions,  it  follows  that  the  velocity  field  due  to  any 
doublet  distribution  may  be  replaced  by  that  due  either  to  the  source 
distributions  of  Theorem  2  or  to  the  vorticity  distributions  of  Theorem  1 
(plus  the  local  term p.  for  points  lying  inside  the  influencing  volume). 

By  rearranging  the  various  effects,  the  following  theorem  can  be  obtained: 

Theorem  3  "The  velocity  field  due  to  a  combined  volume  distribution 

~  of  sources  of  density  ( ot  =  div Ju  )  and  vorticity  of 

density  (  Yv  =  curl  ft- )  throughout  a  volume  Jl  is  identical 
to  that  due  to  a  combined  distribution 3  on  the  surface 
of  that  volume ,  of  surface  sources  of  density  (  oj  =  fat.  -  *i ) 
and  surface  vorticity  of  density  (  Ts  -  -fat  x-n  J,  together 
with  a  local  velocity  increment  equal  to  the  local  value  of  Jx 
for  points  which  lie  inside  the  volume  fL". 


3.3  Normal  Doublet  Distribution  on  a  Closed  Surface 


Referring  again  to  Fig.  la,  consider  the  velocity  induced  at  a  point  P  due 
to  a  continuous  surface  distribution  of  doublets,  of  density  JUL%  ,  normal  to  the 
closed  surface  S. 

Using  (6)  and  (1)  the  potential  induced  at  P  is  given  by: 

(17)  471  4>P  =  ~  ///**  $radQ.(ty-  & 

3 

which  by  virtue  of  Gauss'  theorem  (k) t  becomes: 

(18)  4 n  =  -JJJ*  (f  dH 

JX 

where  f  is  interpreted  as  a  scalar  function  which  varies  continuously 
throughout  A  ,  but  which  has  the  scalar  value  JULS  at  the  surface  S. 

Now  using  (2) : 

Iff  div*  Cf  =  !![  f)  +  Iff  /  Vq  (i) ASl 

By  an  argument  similar  to  that  used  in  3.1.2,  the  last  term  is  equal  to 
zero  for  points  P  which  lie  outside  SI  ,  and  equal  to  the  local  value  -  bTlfi 
for  points  within  SI  . 

Thus  for  points  outside  SI  ,  (18)  becomes: 

(19)  4^  #p  =  (^4/)  cUi 

Jl 

whilst  for  points  within  SI  it  becomes: 

(19a)  A-Tt$?  =  -  jjff  &))  •  f)  dil  +  4TT/P  . 

By  comparing  with  equation  (6)  it  can  be  seen  that  the  volume  integral  in 
equations  (19)  and  (19a)  is  equal  to  the  potential  induced  at  P  by  a 
volume  doublet  distribution  fa  ,  given  by  fa  -  grade  f  ;  this  doublet 
distribution  is  irrotational  since  curl^grad^f  =  0.  We  thus  obtain  the 
following  theorem: 
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Theorem  4  "The  potential  and  velocity  induced  at  any  point  by  an 

irrotational  volume  distribution  of  doublets  of  arbitrary 
density  =  grad  4  (where  4  is  an  arbitrary  scalar 

function  in  that  volume)  is  identical  to  that  induced  by  a 
surface  distribution  of_normal  doublets  on  the  surface  of 
that  volume >  given  by  JU-s  -  -ftrL  (where  denotes  the  value  of 
4  infinitesimally  inside  S  and  ri  denotes  the  unit  outward 
normal  to  that  surface) ,  augmented  when  P  lies  inside  that 
volume j  by  the  local  values  of  and  of  -grad  f  , 
respectively" . 

By  comparing  the  above  Theorem  with  Theorem  la  for  the  case  where  curl^  is 
zero  (as  above)  the  following  Theorem  is  obtained: 

Theorem  5  "The  velocity  induced  at  any  point  by  a  continuous  surface 
distribution  of  normal  doublets /AsrL  on  any  closed  surface 
is  identical  to  that  induced  by  a  surface  vorticity 
distribution  Ys  on  that  same  surface  given  by  rs  =Cprad2/u,)«  ■n 
where  n  denotes  the  local  unit  normal  to  the  surface  and 
the  operator  grad a  indicates  that  only  surface  derivatives 
are  required  (the  normal  component  of  grad^u  contributes 
nothing  to  the  vector  product  with  -n  ) 


3-1*  General  Statement  of  Equivalent  Distributions 

By  comparing  Theorems  2,  h  and  5  in  the  special  case  where  the  volume 
doublet  distribution  /£>  throughout  the  volume  SI  is  such  that  div^yOv  =  O  and 
curla  yJv  =  0,  and  is  defined  by  Jdv-  grad  f,  the  following  general  theorem  is 
obtained: 

Theorem  6  "The  following  distributions  produce  identical  velocity  fields: 

(i)  The  arbitrary  doublet  distribution  JIy=grad  4  such  that  Vlf  =  0 
throughout  Jl  ; 

(ii)  The  surface  source  distribution  <ys  =  -ft .  grad  4  on  the  closed 

surface  S  offl ,  grad  f  being  evaluated  infinitesimally  inside  S; 

(Hi)  The  surface  normal  doublet  distribution  /us=4'non  the  closed  surface 
S  of  Jl j  4  being  evaluated  infinitesimally  inside  S ;  for  points 
lying  inside  SI  the  velocity  due  to  this p-s  must  be  augmented  by  the 
local  value  of  -  grad  4  i 

(iv)  The  surface  vorticity  distribution  3^  =  grad  f  > l  ,  grad  f 

being  evaluated  infinitesimally  inside  S ;  for  points  lying  inside  Si 

the  velocity  due  to  this  Ts  must  be  augmented  by  the  local  value  of 
-grad  f 

Note  that  the  alternative  surface  distributions  may  lie  on  different  surfaces 
of  the  volume.  For  example  volume  doublicity  parallel  to  the  axis  of  a  prism 
may  be  replaced  by  surface  source  distributions  on  the  end-faces  normal  to 
that  axis,  or  by  surface  vorticity  distributions  on  the  remaining  surfaces. 
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4.  APPLICATION  TO  CERTAIN  SIMPLE  CASES 


Some  results  will  now  be  derived  which  are  of  relevance  to  surface  singularity- 
methods. 

4.1  Doublet  Distribution  on  a  non-closed  Surface 

Consider  an  arbitrary  surface  S  with  continuously  varying  tangent  plane.  S 
has  the  perimeter  L.  Consider  the  volume  J1  swept  out  by  displacing  every 
point  of  S  by  a  constant  distance  h  in  the  direction  of  the  normal  to  the 
local  tangent  plane  of  S.  It  is  assumed  that  h  is  less  than  the  smallest  radius 
of  concave  curvature  of  S  so  that  normals  to  S  do  not  intersect  within  fl . 

The  boundary  of  11  can  then  be  defined  by  the  displaced  surface  S'  (which  also 
has  continuously  varying  tangent  plane),  by  the  'closed'  edge  E,  and  by  the 


Consider  now  the  surface  £  passing  through  all  points  Qf  which  are  contained 
within  this  volume  at  a  constant  normal  distance  fh  from  S  (0  <  f  <  1); 
suppose  that  is  defined  by  the  equation  F#(x,y,z)  =  0  where  (x,y,z)  are 
the  cartesian  coordinates  of  a  point  on  this  surface,  and  that  E  is  defined 
by  G(x,y,z)  =  0.  At  any  point  the  vector  N  =  f~aci  Ft  (x,y,z)  defines  a 
vector  which_is  locally  normal  to  £ (and  parallel  to  the  unit  vector  ri  at  the 
point  where  N  intersects  S).  In  general  N  will  not  be  a  unit  vector;  its_ 
magnitude  will  vary  over  Sand  indeed  for  points  lying  along  the  same  line  N' 
but  with  different  values  of  f. 

Suppose  no  that  the  above  volume  Jl  contains  a  doublet  distribution  which  at 
any  point  xS  defined  by  where  m  is  a  scalar  function  which  is  constant 

along  any  normal  to  the  surface  S  but  is  otherwise  variable  throughout  the 
vc 1 ume  fl  .  At  some  external  point  P  which  is  sufficiently  far  removed  fromjl  , 
this  volume  doublet  distribution  will  produce  approximately  the  same  potential  ^ 
and  velocity  as  a  surface  doublet  distribution  of  variable  surface  density  TnKNn 
on  some  'mean'  surface  between  S  and  S',  where  N  is  the  mean  magnitude  of  N 
along  the  local  normal  to  S,  and  n  is  the  local  unit  normal  to  S. 

Now,  according  to  Theorem  1,  this  volume  doublet  distribution  /£,  will  produce 
exactly  the  same  velocity  at  any  external  point  P  as  a  volume  vorticity 
distribution  of  density  %  =  curls/itv  throughout  SI  ,  together  with  a  surface 
vorticity  of  density  %  =  on  the  boundaries  S,  S'  and  E,  n  being  the 

local  outward  unit  normal  for  each  of  these  surface.  Since  is  defined  to 
be  normal  to  S  and  S',  the  equivalent  7^  on  these  surfaces  is  zero.  The 
equivalent  volume  vorticity  distribution  Yy  can  be  written,  using  equation 
(2a):  _  _ 

Ty  =  cuurlq  jlv  =  curln  m  N  =  (grad^  m.)  x  N  +■  m  cuv-C*  N 

-  (yrad.^  m")  x  N 

since  curlft  N  =  curlft  ffradq  2)  =  0. 

It  may  be  noted  that  grad^m  will  have  no  component  in  the  direction  of  ¥ 
since  m  is  defined  to  be  constant  in  that  direction.  Also,  since  h  is 
constant,  we  can  write  (grad^  m)  =  -^(grad,  mh.)  and  this  vector  will  be 
constant  for  all  points  which  lie  on  a  particular  normal  to  S.  Thus  at  an 
external  point  sufficiently  far  removed  from  fl  ,  the  volume  vorticity 
distribution  will  produce  approximately  the  same  velocity  as  a  surface 
vorticity  distribution  of  variable  surface  density  fs  =  (grad2  mhN)  x  n 
on  some  'mean'  surface;  the  operator  gradx  excludes  derivatives  normal  to 
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this  surface;  N  is  again  the  mean  magnitude  of  N  along  the  local  normal 
to  S,  and  n  the  local  unit  normal  to  S. 

The  surface  vorticity  distribution  Ts  on  the  'edge'  surface  E  can  be 
written: 

X  =  /Tv  *  £  =  (*»€  Nft)  X  4*  =  ■£  Njn 

where  n4is  the  local  normal  to  E  and  me  and  Ne  are  the  local  values  of  M  and  N 
at  the  surface  E. 

At  an  external,  point  sufficiently  far  removed  from  JL  this  vorticity 
distribution  on  the  edge  surface  E  produces  approximately  the  same  velocity  A 
as  a  line  vortex  distribution  of  variable  line  density  V  defined  by  f  =  ?n4h  N  t 
where  t  =  n  x  n4  is  the  unit  vector  along  the  perimeter  of  the  'mean'  surface. 

In  the  limit  as  h  decreases  to  zero,  whilst  the  function  m.  is  increased  such 
that  the  function  mh  5?  remains  finite  { jX^  say),  the  'mean'  surface  becomes 
coincident  with  S  and  N  becomes  the  vector  grad  Fc(x,  uln')  pertaining  to  S. 

Consequently  the  above  'approximate'  statements  reduce  to  the  following  exact 
statement : 

"The  velocity  field  induced  by  a  distribution  of  normal  doublets  of  surface 
density  ,  on  any  bounded  surface  S  with  continuously  turning  tangent 
plane,  is  identical  to  that  induced  by  a  surface  vorticity  distribution 
Ts  on  that  same  surface,  defined  by  ss  =  grad 2  JUS  (where  ri  is  the  local 
unit  normal  to  S ),  together  with  a  line  vortex  f  along  the  perimeter  L  of  S  , 
defined  by  f  =  t  ,  where  JU^is  the  magnitude  of  at  that  perimeter, 
and  £  is  the  unit  vector  along  L  such  that  the  vector  £  k  JHS  points  away 
from  S 

This  statement  forms  an  extension  to  Theorem  5»  which  was  derived  for  a  closed 
surface  only. 

Note  that  the  above  statement  has  been  derived  for  any  general  surface  S  with 
continuously  turning  tangent  plane;  however,  it  can  be  seen  that  it  is  also 
valid  for  the  case  of  surfaces  containing  lines  at  which  the  tangent  plane  is 
discontinuous,  provided  that  the  magnitude  of  the  doublet  density  is 
continuous  across  such  a  line;  similarly  the  perimeter  L  need  not  form  a 
line  with  continuously  turning  tangent. 

In  the  special  case  of  surface  doublicity  normal  to  S  and  of  uniform  density  fX. 
over  S,  the  equivalent  surface  vorticity  density  is  zero  (since  grad2/Mf  =  0) 
so  that  this  doublet  distribution  may  be  replaced  by  the  concentrated  edge  vortex 
P  = jut  alone,  this  line  vortex  being  of  constant  line  density.  Conversely, 
it  can  be  seen  that  a  constant-strength  line  vortex  along  any  arbitrary  closed 
contour  produces  a  velocity  field  identical  to  that  produced  by  a  uniform 
normal  doublet  distribution  on  any  arbitrary  surface  having  that  contour  as 
perimeter. 

It  is  also  worth  noting  that  a  surface  containing  doublicity  whose  axis  is 
everywhere  tangential  to  that  surface  and  of  density fx  is  equivalent,  by  a 
similar  argument,  to  a  surface  source,  distribution  of  density  -  divayu 
together  with  a  line  source  of  density  nj  =  p~  ?  where  s  denotes  the  local 
unit  normal  to  the  perimeter  of  S,  lying  in  the  tangent  plane  to  S  at  that 
perimeter  and  orientated  away  from  S. 
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In  the  special  case  of  surface  doublicity  having  its  axis  everywhere  tangential  to S 
and  unidirectional  and  its  density  yx  uniform  over  S,  the  equivalent 
surface  source  density  is  zero,  so  that  this  doublet  distribution  may  be 
replaced  by  a  concentrated  edge  source  of  variable  line  density  =  yU.  s. 
Advantage  will  be  taken  of  this  fact  in  section  4.2. 


4.2  Velocity  Induced  by  a  Uniform  Source  Distribution  on  a  Planar  Polygon 

Consider  a  planar  polygon  S  carrying  a  tangential  surface  doublet 
distribution yU  which  is  unidirectional  and  which  has  the  uniform  density  K. 

Suppose  that  a  cartesian  coordinate  system  (x,y,z)  is  defined  such  that  the 
polygon  lies  in  z  =  0  and  that  the  x  axis  is  parallel  to  the  axis  of  the 
doublicity.  Thus  we  can  writ e yU  ~  K  t .  As  indicated  in  section  4.1,  this 
doublet  distribution  produces  the  same  potential  and  velocity  field  as  a 
concentrated  line  source  of  density  (Ks.i)  along  the  edges  of  the  polygon. 

On  each  straight  edge  of  the  polygon,  the  vector  §  is  constant,  so  that  the 
equivalent  line  source  along  that  edge  is  uniform  in  density. 

It  can  be  seen  by  comparing  equations  (5)  and  (6)  that  the  potent ial  induced 
at  any  point  by  the  doublet  distribution  Ki  on  S  is  the  same  as  the  negative 
of  the  x-component  of  the  velocity  induced  at  that  point  by  a  source 
distribution  of  uniform  density  K  on  S. 

It  follows  that  the  x-component  of  the  velocity  induced  at  any  point  by  a 
uniform  surface  source  distribution  of  density  K  on  S  is  identical  to  the  sum  of 
the  potentials  induced  at  that  point  by  a  line  source  along  each  edge  of  the 
polygon  of  density  -Ks.l  (constant  al^ng  each  edge).  Similarly  the  y-component 
can  be  obtained  from  the  sum  of  the  potentials  due  to  line  sources  of  density  -ltf.j. 
These  potentials  are  obtained  very  simply  by  integrating  the  elemental  expression 
given  in  equation  (5). 

The  z-component  of  the  velocity  may  be  derived  by  considering  the  potential 
due  to  a  uniform  c' sublet  distribution  defined  by  /t  =  fdc  (i.e.  doublets  with 
axes  normal,  to  the  plane  S).  In  this  case  the  equivalent  system  consists  of 
a  line  vortex  of  constant  density  K  on  each  edge  of  the  polygon;  the 
potential  at  P  due  to  this  closed  ring  vortex  may  be  obtained  by  integrating 
(from  infinity  to  the  point  P,  with  respect  to  z)  the  z-component  of  the 
velocity  induced  by  this  vortex.  This  latter  function  is  obtained  by 
integrating  equation  (7),  (or  alternatively  from  the  velocity  induced  by  a 
line  sourcje  of  density  K,  using  the  equivalence  expressed  after  equation  (7)). 

It  can  thus  be  seen  that  the  velocity  induced  by  a  uniform  surface  source 
distribution  on  a  plane  polygon  may  be  obtained  purely  from  line  integrals 
evaluated  along  the  edges  of  that  polygon;  this  fact  is  utilised  in 
'panel  methods'  which  employ  such  source  distributions.  Similar  arguments 
may  be  used  to  simplify  the  integrals  required  in  higher-order  methods 
employing  non-uniform  source  distributions  or  non-planar  panels. 
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4 . 3  Source  distribution  on  the  Surface  of  a  Wedge 


Consider  the  volume  formed  by  rotating  the  rectangle  ABCD  through  an  angle  o( 
about  the  z  axis  which  is  parallel  to  AB  and  a  distance  r.  from  it.  Points 
within  this  wedge  will  be  defined  in  terms  of  the  cylindrical  polar  coord¬ 
inates  (r ,  9  ,2) 

Suppose  that  the  volume  contains  a  doublet  distribution  /I  of  volume  density  jjjt 
where  ja„  is  constant,  the  axis  of  each  elemental  doublet  being  normal  to  the 
plane  containing  itself  and  Oz  (with  the  positive  end  uppermost).  It  can 
easily  be  shown  that  both  curl^Z  and  div JX  are  zero  throughout  the  volume. 

Using  Theorem  1,  the  volume  doublet  distribution  produces  the  same  external 
velocity  as  a  surface  jvorticity  distribution  on  the  surfaces  of  that 
volume,  described  by  %  =  ju*n  (the  volume  vorticity  has  zero  density  since 
curl Ju.  =  0).  Since  the  doublicity  is  locally  normal  to  the  rectangular 
faces,  the  value  of  fs  on  those  faces  is  zero,  and  we  are  left  with  vorticity 
which  has  surface  density  Jju/r  on  the  remaining  faces. 

On  the  surface  St  this  vorticity  vector  points  in  the  direction  of  the  z  axis 
whereas  on  SA  it  points  in  the  opposite  direction,  on  SA  it  points  radially 
away  from  the  z  axis,  and  on  Sg  radially  towards  it.  In  the  limit  as  r  -*  0, 
the  surface  vorticity  distribution  on  S4  ,  of  density  /i./r.  ,  distributed 

over  an  axe  length  of  o(Tt  ,  produces  the  same  velocity  field  as  a  line  vortex 

of  line  density  o(ju.0  ,  along  the  line  AB  which  now  lies  on  the  z  axis. 

Using  Theorem  2,  on  the  other  hand,  the  volume  doublet  distribution  produces 
the  same  external  velocity  as  a  surface  source  distribution  on  the  surfaces 
of  that  volume,  of  surface  density  Oj  =  jL  .n  (the  volume  source  density  is 
zero  since  div  jui  =  0).  This  surface  source  density  is  zero  except  on  the 
two  rectangular  plane  faces,  on  which  it  has  the  variable  magnitude  . 

The  source  density  is  positive  on  the  upper  surface  and  negative  on  the  lower 
surface. 

Considering  the  case  where  ABCD  is  of  indefinite  extent  and  r0  =  0,  so  that 

the  volume  corresponds  to  an  infinite  wedge  of  included  angle  o< ,  it  can  be 

seen  that  the  velocity  field  induced  outside  the  wedge  by  a  line  vortex  of 
density  T  along  its  ridge  corresponds  to  that  induced  by  surface  source  and 


sink  distributions  on  its  two  faces,  of  surface  density  <f  =  T  /^f  ;  this 
density  clearly  approaches  infinity  near  the  ridge  (i.e.  as  r  -*0). 

By  adopting  the  results  derived  in  section  4.1,  the  line  vortex  may  be 
replaced  by  a  uniform  surface  distribution  of  normal  doublets  over  any 
indefinitely  large  surface  which  terminates  at  the  ridge  of  the  wedge. 

These  results  (or  modified  versions  of  them,  using  line  vortices  which  are 
not  of  constant  strength,  or  surfaces  more  complex  than  an  infinite  wedge) 
are  of  relevance  in  problems  where  wakes,  modelled  by  doublet  sheets,  arise 
from  wing  trailing  edges  or  from  smooth  surfaces  (wedge  angle  o(  =  TT  )  which 
are  modelled  by  surface  source  distributions. 


4.4  Modelling  of  Flow  Past  an  Infinite  Circular  Cylinder;  Jets 

tz 


solid  cylinder  (a)  2-D  Doublet 

It  is  well  known  that  the  potential,  flow  outside  a  solid  2-D  circular  __ 
cylinder  of  radius  a  installed  perpendicular  to  a  uniform  free  stream  U„ 
can  be  modelled  by  the  interaction  between  li„  and  a  2-D  line  doublet  of_ 
line  density  27Ta2H„<(  ,_the  doublet  axis  being  opposite  in  direction  to  UM 
(as  sketched  above);  U*,  is  assumed  parallel. to  the  x  axis  with  x  =  0  at  the 
centre  of  the  cylinder. 


It  is  also  easy  to  demonstrate  that  the  cylinder  may  be  replaced  by  a  surface 
source  distribution  of  surface  density  oj  =-211,0  cos  0  : 


By  Theorem  6  it  can  then  be  seen  that  this  source  distribution  will  induce  the 
same  external  velocity  as  a  volume  doublet  distribution  given  by 
fl-v  =  grad  (-2U„  x)  =  -  2U^  1,  i.e.  a  uniform  volume  distribution  of  doublets 
with  the  doublet  axis  pointing  in  the  -x  direction: 
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This  volume  doublet  distribution  is  also  equivalent  to  either  a  surface 
distribution  of  vorticity  =  2U^  sin  $  or  of  normal  doublets  of  surface 
density  J-<s=  -2U^  a  cos  9  : 


In  this  special  case,  the  volume  doublet  distribution  (c)  is  equivalent  to 
either  any  one  of  the  surface  distributions  (b,  d,  e)  or  the  line 
distribution  (a). 


Similar  arguments  may  be  extended  to  cases  in  which  the  cylinder  is  not  'solid' 
or  'straight'  or  'circular';  such  cases  arise  in  the  modelling  of  jets 
issuing  obliquely  into  a  uniform  stream.  Such  jets  take  up  a  curved  shape 
along  their  axis,  and  have  a  non-circular  cross-section;  fluid  is  entrained 
into  the  jet  through  its  boundary;  the  inside  of  the  jet  contains  a  pair  of 
contra-rotating  vortices. 

The  optimum  modelling  of  this  jet  is  a  current  subject  of  research  and  will 
not  be  considered  in  detail  here,  but  it  is  worth  noting  that  the  equivalent 
representations  of  Theorem  6  may  allow  certain  of  the  models  under  development 
to  be  considerably  reduced  in  terms  of  their  computational  requirement. 
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CONCLUSIONS 


Certain  theorems  have  been  derived  which  state  the  equivalence  of  particular 
line,  surface  and  volume  distributions  of  sources,  doublets  and  vorticity. 
These  theorems  allow  the  influence  of  singularity  distributions  employed,  for 
example^  in  'panel  methods',  to  be  expressed  in  terms  of  the  influence  of  more 
elementary  distributions.  It  is  possible  that  these  results  will  be  of  some 
value  in  the  computation  of  compressible  and/or  viscous  flows  and  in  the 
modelling  of  non-linear  problems  involving,  for  example,  thick  wakes  or 
lifting  jets. 
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